Acta Polytechnica Vol. YY No. X/ZZZZ 



ARGO-YBJ: Status and Highlights 



Giuseppe Di Sciascio^ on behalf of the ARGO-YBJ Collaboration 



(N 

O 

(N 
> 

o 

O: 

^' 

6" 
^ . 

a. 



> 
m 

d 



X: 



^ INFN - Sezione di Roma Tor Vergata, Viale della Rtcerca Scientifica 1, 1-00133 Roma 
Corresponding author: disciascio@roma2.infn.it 

Abstract 

The ARGO-YBJ experiment is in stable data taking since November 2007 at the YangBaJing Cosmic Ray Laboratory 
(Tibet, P.R. China, 4300 m a.s. 1., 606 g/cm^). ARGO-YBJ is facing open problems in Cosmic Ray (CR) physics in 
different ways. The search for CR sources is carried out by the observation of TeV gamma-ray sources both galactic and 
extra-galactic. The CR spectrum, composition and anisotropy are measured in a wide energy range (TeV - PeV) thus 
overlapping for the first time direct measurements. In this paper we summarize the current status of the experiment and 
describe some of the scientific highlights since 2007. 
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1 Introduction 

Exploiting the full coverage approach at high alti- 
tude, the ARGO-YBJ experiment is an air shower 
array able to detect the cosmic radiation with an en- 
ergy threshold of a few hundred GeV. The detector 
is in stable data taking since November 2007 with a 
duty cycle larger than 86%. The trigger rate is 3.5 
kHz. The detector characteristics and performance 
are described in detail in [U [2]. The main results 
obtained by ARGO-YBJ are described in [7. In this 
paper the status of the experiment is summarized and 
some highligths reviewed. 



2 Gamma- ray astronomy 

From November 2007 the ARGO-YBJ experiment 
collected about 4x10^^ events in 1543 days of total 
effective observation time. Five known VHE 7-ray 
sources have been detected with a statistical signif- 
icance greater than 5 standard deviations (s.d.), i.e. 
Crab Nebula, Mrk421, Mrk501, MGRO J1908-F06 
and MGRO J203H-41. A number of flares from 
Mrk421 and Mrk501 were observed and studied in 
detail. Evidence of a TeV flaring activity from the 
Crab Nebula in coincidence with AGILE/Fermi ob- 
servations is also reported. Details on the analysis 
procedure (e.g., data selection, background evalua- 
tion, systematic errors) are discussed in [H [S] . 

2.1 Crab Nebula 

With ah the data recorded in 3.5 years ARGO-YBJ 
observed a TeV signal from the Crab Nebula with a 
statistical significance of about 17 s.d., proving that 



the cumulative sensitivity of the detector reached 0.3 
Crab units. 

The observed flux is consistent with a steady 
emission, and the observed differential energy spec- 
trum in the 0.3-30 TeV range can be described 
by dN/dE = (3.0±0.30) x 10"" (E/TeV)-2 "±o.o9 
photons cm^^ s^^ TeV^^, in good agreement with 
other observations. According to MC simulations, 
84% of the detected events comes from primary pho- 
tons with energies greater than 300 GeV, while only 
8% comes from primaries above 10 TeV. We evaluate 
a systematic error on the flux less than 30% mainly 
due to the background estimaation and to the uncer- 
tainty on the absolute energy scale. 

According to the AGILE and Fermi data 4 ma- 
jor flaring episodes at energies E >100 MeV occurred 
during the ARGO-YBJ data acquisition [i [3 [HI [H] • 

Flare 1: starting time MJD 54857 (Feb. 2009), 
duration At'^16 days, maximum flux Ymax ^5 times 
larger than the steady flux f^- During this flare no 
excess is present in our data, for any multiplicity 
threshold. 

Flare 2: starting time MJD 55457 (Sept. 2010), 
duration At ^4 days, maximum flux ¥max '^5 times 
larger than the steady flux [SI [7] . According to tem- 
poral data analysis the 7-ray emission is concentrated 
in 3 narrow peaks of ^12 hours duration each [l0l[8]. 
Integrating the 3 transits we observe for Npad >40 
an excess of 3.1 s.d. over the expected steady flux 
(0.55 s.d.). If the excess were due to a flare, the 
7-ray flux would be higher by a factor ~5 with re- 
spect to the steady flux at energies around 1 TeV. 
Integrating the data over 10 transits (from MJD 
55456/57 to MJD 55465/66) the signal signiflcance is 
4.1 s.d. (pre-trial), while 1.0 s.d. is expected from the 
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steady flux [IT. No measurement from Cherenkov 
telescopes is available in coincidence with our ob- 
servations and with the different spikes to confirm 
this excess. Sporadic measurements performed by 
MAGIC and VERITAS telescopes at different times 
from MJD 55456.45 to MJD 55459.49 show no evi- 
dence for a flux variability [HI [T^ . 

Flare 3: starting time MJD 55660 (Apr. 2011) 
[14| . duration At '^6 days, maximum flux F^ax ^14 
times larger than the steady flux f8| . Integrating the 
ARGO-YBJ data over the 6 days in which AGILE de- 
tected a flux enhancement, i.e. from MJD 55662 to 
MJD 55668, we find evidence for an excess for events 
with Npad > 100 at a level of about 3.5 s.d. . No 
measurement from Cherenkov telescopes is available 
during these days, due to the presence of the Moon 
during the Crab transits. 

Flare 4: starting time MJD 56111 (July 3rd, 
2012) duration At '--^3 days. The daily-averaged 
emission doubled from (2.4±0.5)x 10~^ ph cm~^ 
sec~^ on July 2nd to (5.5±0.7) x 10—6 ph cm^^ sec^^ 
on July 3rd. A preliminary analysis of the ARGO- 
YBJ data shows an excess of events with statistical 
significance of about 4 s.d. from a direction consis- 
tent with the Crab Nebula on July 3rd, corresponding 
to a flux « 8 times higher than the average emission 
at a median energy of ~1 TeV [H]. The expected 
steady flux corresponds to 0.33 s.d. . No significant 
exces is detected in the following days from July 4th 
to 6th. Once again, no measurement from Cherenkov 
telescopes is available during these days. 

In conclusion, the ARGO-YBJ data show a 
marginal evidence of a TeV flux increase correlated to 
the MeV-GeV flaring activity, with a statistical sig- 
nificance not enough to draw a firm conclusion (the 
post-trial chance probability is of order 10~^ for each 
flare). Nevertheless, the chance probability of observ- 
ing 3 flares out of 4 in coincidence with the satellites 
is greatly depressed. A detailed analysis of the Crab 
Nebula TeV emission is under way and a paper is in 
preparation. 



2.2 The blazar Mrk421 

ARGO-YBJ is monitoring Mrk421 above 0.3 TeV, 
studying the correlation of the TeV flux with X-ray 
data. We observed this source with a total signifi- 
cance of about 14 s.d., averaging over quiet and ac- 
tive periods. As it is well known, this AGN is char- 
acterized by a strong flaring activity both in X-rays 
and in TeV 7-rays. Many flares are observed simul- 
taneously in both bands. The 7-ray flux observed by 
ARGO-YBJ has a clear correlation with the X-ray 
flux. No lag between X-ray and 7-ray photons longer 
than 1 day is found. The evolution of the spectral en- 
ergy distribution (SED) is investigated by measuring 
spectral indices at four different flux levels. Spec- 



tral hardening is observed in both X-ray and 7-ray 
bands. The 7-ray flux increases quadratically with 
the simultaneously measured X-ray flux. All these 
observational results strongly favor the Self-Synchro 
Compton (SSC) process as the underlying radiative 
mechanism. The results of Mrk421 long-term moni- 
toring are summarized in |16j . 

2.3 The blazar MrkSOl 

The long-term observation of the Mrk501 TeV emis- 
sion by ARGO-YBJ can be described by the following 
differential energy spectrum: dN/dE = (1.92±0.44) 
X 10-12 (E/2 TeV)-2-59±°-27 photons cm-^ 

TeV~\ corresponding to 0.312±0.076 Crab units 
above 1 TeV. 

In October 2011 started the largest flare since 
2005. During the brightest 7-ray flaring episodes 
from October 17 to November 22, 2011, an ex- 
cess of the event rate over 6 s.d. has been de- 
tected by ARGO-YBJ, corresponding to an increase 
of the 7-ray flux above 1 TeV by a factor of 6.6±2.2 
from its steady emission [T7|. During the flare, the 
differential flux is dN/dE = (2.92±0.52) x lO'^^ 



(E/4 TeV) 



-2.07±0.21 



photons cm s TeV , cor- 



responding to 2.05±0.48 Crab units above 1 TeV. 
Remarkably, 7-rays with energies above 8 TeV are 
detected with a statistical significance of about 4 s.d., 
which did not happen since the 1997 flare. The aver- 
age SED for the steady emission is well described by a 
simple one-zone SSC model. However, the detection 
of gamma-rays above 8 TeV during the flare chal- 
lenges this model due to the hardness of the spectra 
[17]. 
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Figure 1: Gamma-ray flux from MGRO J1908+06 
measured by ARGO-YBJ (red line) compared to 
other measurements. The dashed area represents the 
1 s.d. error. The plotted errors are purely statistical 
for all the detectors. For details and references see 
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2.4 MGRO J1908+06 

The 7-ray source MGRO J1908+06 was discovered 
by Milagro at a median energy of ■~20 TeV and con- 
firmed by HESS at energies above 300 GeV. The Mi- 
lagro and HESS energy spectra are in disagreement, 
being the Milagro result about a factor 3 higher at 
10 TeV. 

ARGO-YBJ observed a TeV emission from 
MGRO J1908+06 with a maximum significance of 
about 7.3 s.d. for Npad >20 in 6867 hours on-source 
[18] ■ The intrinsic extension is determined to be 
a^xt = 0.49°±0.22, consistent with the HESS mea- 
surement {<7ext — 0-34°lg Q3). The best fit spec- 
trum obtained is: dN/dE = (6.1±1.4) x 10"" (E/4 
rpgY)-2.54±o.36 photons cm'^ s-i TeV-\ in the en- 
ergy range 1-20 TeV (see Fig. [T|). The measured 
gamma-ray fiux is consistent with the Milagro re- 
sults, but is ~2-3 times larger than the flux derived 
by HESS at energies of a few TeV. Given the reduced 
significance of the excess at high energies, we are not 
able to constrain the shape of the spectrum above 
10 TeV and to definitively rule out a possibile high 
energy cutoff. 

The continuity of the Milagro and ARGO-YBJ 
observations and the stable excess rate observed 
by ARGO-YBJ along 4 years of data taking sup- 
port the identification of MGRO J1908-I-06 as a sta- 
ble extended source, likely the TeV nebula of PSR 
J1907-h0602, with a flux at 1 TeV -67% that of the 
Crab Nebula. Assuming a distance of 3.2 kpc, the in- 
tegrated luminosity above 1 TeV is ~1.8 times that 
of the Crab Nebula, making MGRO J1908-(-06 one 
of the most luminous Galactic gamma-ray sources at 
TeV energies [18]. 



sion. The source VER J2016-I-372 is marked with a 
small circle without position errors. The small circle 
within the errors of MGROJ203H-41 indicates posi- 
tion and extension of the source TeV J2032-I-4130 as 
estimated by the MAGIC collaboration. The open 
stars mark the location of the 24 GeV sources in the 
second Fermi LAT catalog. For details and refer- 
ences see [in]. 



2.5 MGRO J2031+41 and the 
Cygnus region 

The Cygnus region contains a large column density 
of interstellar gas and is rich of potential CR accel- 
eration sites as Wolf-Rayet stars, OB associations 
and supernova remnants. Several VHE gamma-ray 
sources have been discovered within this region in the 
past decade, including two bright extended sources 
detected by the Milagro experiment. 
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Figure 2: Signiflcance map of the Cygnus region as 
observed by ARGO-YBJ for Npad >20. The four 
known VHE 7-ray sources are reported. The er- 
rors on the MGRO source positions are marked with 
crosses, while the circles indicate their intrinsic sizes. 
The cross for VER J2019-h407 indicates its exten- 



Figure 3: Energy spectrum of MGRO J2031+41 /TeV 
J2032-h4130 as measured by ARGO-YBJ (magenta 
solid line). The spectral measurements of HEGRA 
and MAGIC are also reported for comparison. The 
black solid line and shaded area indicate the differ- 
ential energy spectrum and the 1 s.d. error region as 
recently determined by the Milagro experiment. The 
two triangles give the previous flux measurements by 
Milagro at 20 TeV and 35 TeV. For details and ref- 
erences see [l9] . 

The gamma ray source MGRO J2031-I-41, de- 
tected by Milagro at a median energy of ~ 20 
TeV, is spatially consistent with the source TeV 
J2032-h4130 discovered by the HEGRA collaboration 
and likely associated with the Fermi pulsar IFGL 
J2032. 2+4127. The extension measured by Mfla- 
gro, 3.0°±0.9°, is much larger than that initially es- 
timated by HEGRA (about 0.1°). 

The bright unidentified source MGRO J2019+37 
is the most significant source in the Milagro data set 
apart from the Crab Nebula. This is an enigmatic 



Acta Polytechnica Vol. YY No. X/ZZZZ 



source due to its high flux not being confirmed by 
other VHE gamma-ray detectors. Recently, a deep 
survey carried out by VERITAS with a sensitivity 
of ~1% Crab units showed a complex emitting re- 
gion with different faint sources inside the MGRO 
J2019+37 extension. The estimated flux is much 
weaker than that determined by Milagro. 

The Cygnus region has been studied by ARGO- 
YBJ with data collected in a total effective obser- 
vation time of 1182 days [19]. The results of the 
data analysis are shown in Fig. [2] and Fig. |3l A 
TeV emission from a position consistent with MGRO 
J2031-h41/TeV J2032-f4130 is found with a signifi- 
cance of 6.4 s.d. . Assuming the background spec- 
tral index -2.8, the intrinsic extension is determined 
to be Gext = (0.2;to'2)°i consistent with the estima- 
tion by the MAGIC and HEGRA experiments, i.e., 
(0.083±0.030)° and (0.103±0.025)°, respectively. 

The differential fiux (Fig. [S]) is dN/dE = 
(1.40±0.34) X 10~" (E/TeV)-2-^3±o.37 photons 



TeV" 



^, in the energy range 0.6 - 7 TeV. 
= 0.1 the integral flux is 31% that of 



Assuming a^xt 
the Crab at energies above 1 TeV, which is higher 
than the flux of TeV J2032-h4130 as determined by 
HEGRA (5%) and MAGIC (3%). Again, this mea- 
surement is in fair agreement with the Milagro result. 

The reason for the large discrepancy between the 
fluxes measured by Cherenkov telescopes and EAS 
arrays (ARGO-YBJ and Milagro) is still unclear. 
Possible contributions from diffuse 7-ray emission, 
nearby sources, and systematic uncertainties are not 
enough to explain the discrepancy [19^. 

No evidence of a TeV emission above 3 s.d. is 
found at the location of MGRO J2019-H37 and flux 
upper limits at 90% c.l. are set. At energies above 
5 TeV the ARGO-YBJ exposure is still insufficient 
to reach a firm conclusion while at lower energies the 
ARGO-YBJ upper limit is marginally consistent with 
the spectrum determined by Milagro. The observa- 
tion by ARGO-YBJ is about five years later than 
that by Milagro. A flux variation over the whole ex- 
tended region cannot be completely excluded. If the 
flux variation were dominated by a smaller region in 
the source area, the picture could be more reason- 
able. In such a scenario, however, identifying MGRO 
J2019-I-37 as a PWN could be a dilemma because it 
otherwise should have a steady flux. 



3 Cosmic Ray Physics 

Several interesting results have been obtained by 
ARGO-YBJ in the CR physics, as discussed in [3]. 
In the following sections the measurements of the 
anisotropy in the CR arrival direction distribution 
and of the light component (p+He) spectrum of CRs 
are described. 



3.1 Large Scale CR Anisotropy 

The observation of the CR large scale anisotropy by 
ARGO-YBJ is shown in Fig. |4] as a function of the 
primary energy up to about 25 TeV. 
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Figure 4: Large scale CR anisotropy observed by 
ARGO-YBJ as a function of the energy. The color 
scale gives the relative CR intensity. 
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Figure 5; Amplitude of the first harmonic as a func- 
tion of the energy, compared with other measure- 
ments. 
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Figure 6: Phase of the first harmonic as a function 
of the energy, compared with other measurements. 

The so-called Hail-in' and 'loss-cone' regions, cor- 
related to an enhancement and a deficit of CRs, re- 
spectively, are clearly visible with a statistical signif- 
icance greater than 20 s.d.. The tail- in broad struc- 
ture appears to break up into smaller spots with in- 
creasing energy. In order to quantify the scale of the 
anisotropy we studied the 1-D R.A. projections inte- 
grating the sky maps inside a declination band given 
by the field of view of the detector. Therefore, we 
fitted the R.A. profiles with the first two harmon- 
ics. The resulting amplitude and phase of the first 
harmonic are plotted in Fig. [5] and Fig. [6] where 
are compared to other measurements as a function of 
the energy. The ARGO-YBJ results are in agreement 
those of other experiments, suggesting a decrease of 
the anisotropy first harmonic amplitude at energies 
above 10 TeV. 

3.2 Medium Scale Anisotropy 

Fig. [7] shows the ARGO-YBJ sky map in equato- 
rial coordinates containing about 2-10^^ events re- 
constructed with a zenith angle <50°. 
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Figure 7: Medium scale CR anisotropy observed by 
ARGO-YBJ. The color scale gives the statistical sig- 
nificance of the observation in standard deviations. 



The zenith cut selects the declination region b ~ 
-20°-^ 80°. According to simulations, the median en- 
ergy of the isotropic CR proton flux is Ep° wl.8 TeV 
(mode energy wO.7 TeV). 

The most evident features are observed by 
ARGO-YBJ around the positions a - 120°, b - 40° 
and a ~ 60°, 8 ~ -5°, positionally coincident with 
the excesses detected by Milagro [30]. These regions, 
named "1" and "2", are observed with a statistical 
significance of about 14 s.d.. The deficit regions par- 
allel to the excesses are due to a known effect of the 
analysis, that uses also the excess events to evaluate 
the background, overestimating this latter. 

The left side of the sky map seems to be full of 
few-degree excesses not compatible with random fluc- 
tuations (the statistical significance is more than 6 
s.d. post-trial). The observation of these structures 
is reported here for the first time and together with 
that of regions 1 and 2 it may open the way to an 
interesting study of the TeV CR sky. 

To figure out the energy spectrum of the excesses, 
data have been divided into five independent shower 
multiplicity sets. The number of events collected 
within each region is computed for the event map 
(Ev) as well as for the background one (Bg). The 
relative excess (Ev-Bg) /Bg is computed for each mul- 
tiplicity interval. 
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Figure 8: Size spectrum of the regions 1 and 2. 
The vertical axis represents the relative excess (Ev- 
Bg)/Bg. The upper scale shows the corresponding 
proton median energy. The shadowed areas repre- 
sent the 1(7 error band. 

The result is shown in Fig. |8l Region 1 seems 
to have a spectrum harder than isotropic CRs and 
a cutoff around 600 shower particles (proton median 



energy E 



50 _ 



TeV). On the other hand, the excess 



hosted in region 2 is less intense and seems to have a 
spectrum more similar to that of isotropic CRs. We 
note that, in order to filter the global anisotropy, we 
used a method similar to that used by Milagro and 
Icecube. Further studies using different approaches 
are under way. 
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3.3 Light component (p+He) 
spectrum of CRs 
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Figure 9: Light component (p+Hc) spectrum of pri- 
mary CRs measured by ARGO-YBJ compared with 
other experimental results. 

Requiring quasi- vertical showers {9 <30°) and apply- 
ing a selection criterion based on the particle density, 
a sample of events mainly induced by protons and 
helium nuclei with shower core inside a fiducial area 
(with radius '^28 m) has been selected. The con- 
tamination by heavier nuclei is found negligible. An 
unfolding technique based on the Bayesian approach 
has been applied to the strip multiplicity distribution 
in order to obtain the differential energy spectrum of 
the light component (p -I- He nuclei) in the energy 
range (5 - 200) TeV [23]. The spectrum measured 
by ARGO-YBJ is compared with other experimen- 
tal results in Fig. |9l Systematic effects due to dif- 
ferent hadronic models and to the selection criteria 
do not exceed 10%. The ARGO-YBJ data agree re- 
markably well with the values obtained by adding up 
the p and He fluxes measured by CREAM both con- 
cerning the total intensities and the spectral index 
[21] . The value of the spectral index of the power- 
law fit to the ARGO-YBJ data is -2.61±0.04, which 
should be compared with 7p -2.66±0.02 and ^He = 
-2.58±0.02 obtained by CREAM. The present anal- 
ysis does not allow the determination of the individ- 
ual p and He contribution to the measured flux, but 
the ARGO-YBJ data clearly exclude the RUNJOB 
results [22]. We emphasize that for the first time di- 
rect and ground-based measurements overlap for a 
wide energy range thus making possible the cross- 
calibration of the different experimental techniques. 

4 Conclusions 

The ARGO-YBJ detector exploiting the full cover- 
age approach and the high segmentation of the read- 
out is imaging the front of atmospheric showers with 
unprecedented resolution and detail. The digital and 



analog readout will allow a deep study of the CR phe- 
nomenology in the wide TeV - PeV energy range. The 
results obtained in the low energy range (below 100 
TeV) predict an excellent capability to address a wide 
range of important issues in Astroparticle Physics. 

5 Dedication 

This article is dedicated to the memory of my father 
who passed away on September 11th, 2012. 
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